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Abstract 

We analytically study bond percolation on hyperbolic lattices obtained by tiling a hyperbolic 
plane with constant negative Gaussian curvature. The quantity of our main concern is pc2, the value 
of occupation probability where a unique unbounded cluster begins to emerge. By applying the 
substitution method to known bounds of the order-5 pentagonal tiling, we show that Pc2 > 0.382 508 
for the order-5 square tiling, pc2 > 0.472 043 for its dual, and Pc2 > 0.275 768 for the order-5-4 
rhombille tiling. 
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I. INTRODUCTION 



Hyperbolic geometry is an important model of non-Euclidean geometry where mathe- 
maticians have devoted a great deal of efforts since Carl Friedrich Gauss [Ij. In the context 
of statistical physics, hyperbolic geometry has served as a conceptual setting to understand 
geometric frustrations in glassy materials |2|. It also has a nontrivial connection to the 
two-dimensional conformal field theory (s!] and there have been attempts to identify it with 
the underlying geometry of complex networks J]. For this reason, basic understanding of 
physical processes on this geometry is expected to be relevant in a wider context of statis- 
tical physics as well. If we are to discretize a surface by means of regular tiling to study 
physical systems defined on a lattice, in particular, hyperbolic geometry provides infinitely 
more possibilities than the Euclidean geometry: Let {p, q} denote tiling where q regular 
p-gons meet at each vertex. This bracket representation is called the Schlafli symbol. It 
is easy to see that a fiat plane admits only three possibilities: {3,6} (triangular), {4,4} 
(square), and {6, 3} (honeycomb), while every {p, q} such that (p — 2)(g — 2) > 4 describes 
a hyperbolic plane with constant negative Gaussian curvature. In other words, each pair 
of such {p, q} defines a hyperbolic lattice that can completely cover the infinite hyperbolic 
plane with translational symmetry. The most important physical property of a hyperbolic 
plane is that the area of a circle on it is an exponential function of the radius, which means 
that the circumference increases exponentially, too. Therefore, choosing any finite domain 
on a hyperbolic lattice, we find that the vertices at the boundary always occupy a finite 
portion of the whole number of vertices inside the domain even if the domain is very large. 
This property is called nonamenable in literature 5| and makes essential differences in many 
physical systems from their planar counterparts. 

Percolation is a simple yet most interesting problem of ful 
possibility of a global connection out of local connections 
percolation problem, which will be studied in this work: For a given structure of sites and 
bonds linking them, suppose that each bond is open with probability p and closed with 
1 — p, where the parameter p is called occupation probability. One fundamental question 
in percolation is to find a critical value p = Pc where an unbounded cluster of open bonds 
begins to form. This question has been already answered for the three regular ways of tiling 
a fiat plane {t-^ and also for more general ones provided that they allow a generalized cell- 



ly geometric nature, asking the 
6|. Let us introduce the bond 



dual-cell transformation 10|. On these flat lattices, there exists a unique pc above which the 
largest cluster occupies a flnite fraction of the system, and the length scale of this cluster 
becomes unbounded at this point. On a hyperbolic plane, on the other hand, studies of 



percolation started about one decade ago |5 



The most remarkable prediction here is 



that there generally exist two different percolation thresholds Pd and Pc2 with pd < Pc2, 
so that an unbounded cluster begins to appear at pd while a unique unbounded cluster 
is observed only when p reaches a higher value, Pc2- Note that this is a consequence of 
the nonamenable property js] and that these two thresholds coalesce on a flat plane by 
Pel = Pc2 = Pr- Numerical calculations have qualitatively supported this mathematical 



prediction 12|, but a direct numerical estimate of Pc2 is usually a difficult task since the 
system size increases exponentially as the length scale grows. This has led to a debate about 
Pc2 on some hyperbolic structures |l3l-ll5|. 



Recently, nontrivial u ppe r bounds of Pd for self-dual tiling {m, m} were derived by a 
combinatorial argument [16, 17|. If m = 5, for example, Pd'^^ is bounded as 1/4 < Pd'^^ < 



0.381 296. The upper bound is a solution of the following polynomial equation: 

-2 + 3p + 3p2 - 567^=^ + 6721p^ - 35 655/ 
+115 505/ - 257 495/ + 41 8210/ - 509 lOO/ 



+469 900p 
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328 480j9^^ + 171 560p^2 



65 000j9 



13 



+16 900p^^ - 2700p'° + 200p'^ = 0, 
and the lower bound originates from the simple fact that 

Pel > l/(n - 1) 



15 



,16 



(1) 



{5,5} 



< 



for a lattice {m, n} with coordination number n. Gu and Ziff have suggested that Pd 
0.346 and ^^2'^^ > 0.666, with estimating Pd'^^ ~ 0.263 and ^^2'^"^ ~ 0.749 by numerically 
calculating the crossing probability [l^, and these results also strongly support the above 



16|, 



17j is important in two 



analytic bounds. From our point of view, the work in Refs. 
aspects: First, it showed the possibility of rigorous analytic bounds free from any numerical 
ambiguities. Second, it dealt with the problem from a new point of view, that is, in terms 
of the capacity of a quantum erasure channel. In this Brief Report, we point out that 

171] also imply nontrivial bounds of other hyperbolic lattices. 



the nontrivial bounds in Ref. 



Specifically, it is made possible by using the substitution method 



18| . and the lattices 




FIG. 1. Substitution regions of the (a) star-triangle, (b) star-square, and (c) star-pentagon trans- 
formations. The solid lines are occupied with probability p and the dashed lines are occupied with 
Q- 

considered here will be endowed with transitivity to remove any undesirable boundary effects, 
which allows us to exploit duality properties among the lattices, too. In the next section, 
we will briefly explain the substitution method and then show how to apply it to hyperbolic 
lattices as well as the results in Sec. IIIII This work is summarized in Sec. [IVl 



II. SUBSTITUTION METHOD 



The easiest way to explain the substitution method is to begin with the star-triangle 



transformation 



19| . Then, we proceed to other cases such as star-square and star- 



pentagon transformations, which will be used in our problem. Consider a lattice C with 
congruent n-gons as its basic building blocks. By drawing an ra-star with ra-bonds inside 
every n-gon, we obtain another lattice C, where the occupation probability is denoted as q 
in order to avoid confusion with p of C. 



A. Triangle 

If = 3, £ is the triangular lattice [Fig. E^a)], whereas if = 4, £ is the square lattice 
[Fig. m^b)]. To explain the substitution method in a simple manner, we consider the n = 3 
case. Suppose that we happen to know the percolation threshold qc of C. We wish to find 
bounds of Pc on C from the knowledge of qc- Consider one of the triangles T in £ and its 
corresponding star T' in C so that T and T' share the boundary vertices. A, B, and C. We 
look at all the possible cases of connection among the boundary vertices A, B, and C on 
T and T', respectively. Suppose that £ is in the supercritical state and £' is at the critical 
state (percolating phase). Then, we can make the following qualitative statement: it is 
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probable that T has more connectivity among its boundary vertices than T'. To transform 
this quahtative statement into a quantitative expression, we introduce some combinatorial 
concepts. Let iS be a set of all the possible partitions of boundary vertices A, B, and C: If 
A and B are connected by open bonds and C is separated from them, the representation 
is partition AB\C . A set of connected vertices in a partition will be called a block. For 
partition AB\C ^ AB and C are two distinct blocks. And for n boundary vertices in general, 
there are rath order Bell numbers of elements in S [19]. The set S can be a partially ordered 
set if we define an order: For two partitions vr and a of £, we define vr < a if and only if for 
a block 6^ of TT, there exists a block h„ containing 6,^ on a. And we say that vr is more refined 
than cr, or that vr is a refinement of a. For example, < AB\C and AB\C < ABC. 

But there is no order between AB\CD and ABC\D because there is no block in AB\CD 
which covers ABC of ABC\D and vice versa. A subset W of 5 is an up-set if and only if 
for 7ri,7r2 G iS, tti G W and tti < 7r2, then G U. By this definition, we can generate nine 
up-sets as follows: 



Wo 


= S 


Ui 


= {ABC} 


U2 


= {ABC,A\BC} 


W3 


= {ABC,AB\C} 


W4 


= {ABC,B\AC} 


W5 


= {ABC,A\BC,AB\C} 


We 


= {ABC,AB\C,B\AC} 


U-j 


= {ABC,A\BC,B\AC} 




= {ABC, A\BC, AB\C, B\AC} 



Let Pp{U) and Qq{U) be probabilities that T and T' form an element partition of up-set W, 
respectively. Then, we rewrite the qualitative statement as such PpiU) > Qq^iU) for every 
up-set U of S with probability 1. We can solve this inequality with respect to p and get 
a solid interval of Pm < P < 1, for PpiU) is known as a monotone increasing polynomial 
function of p with Pp^^iU) = and Pp=i(U) = 1 for every up-set U ^ S while Qq^iU) is a 
constant. Here, Pm is a lower bound of the percolation threshold of C since if p is smaller 
than p^, the above inequality cannot hold and £ cannot be in the supercritical state. On the 
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contrary, suppose that L is in the subcritical state and CJ is at the critical state. Similarly, 
we can conclude Ppipl) < Qq^iU). By solving this inequality again, we get another interval, 
< p < Pmi where pm is an upper bound of the percolation threshold of C By taking the 
intersection of the two intervals, we obtain < P < Pm- This is the interval in which the 
percolation threshold of C can exist. In n = 3 case, the resulting set of inequalities with 
respect to all its up-sets is found as Pc[Wi] < PcfWj] with i = 0, . . . , 8. In fact, by symmetry 
and triviality, seven of them turn out to be redundant or trivial so we are left with 

3p^(l -p) +p^ < ql, 

3p(l - pY + 3p2(l - p) + p3 < 3^2(^ _ ^ ^3^ 

with < p < 1. The largest value satisfying all these inequalities for given is a lower 
bound of pc- In order to find an upper bound for pc, we need to revert both the inequalities 
above and the smallest p satisfying the reversed inequalities gives us an upper bound of pc- 
The results are shown in Fig. |2]^a). There is one point where the upper and lower bounds 
coalesce, that is, (qcPc) = (l — 2 sin y|, 2 sin y|) • Here, the star-triangle transformation 
yields the exact percolation threshold pc = 2sin(7r/18) for the triangular lattice {3,6} and 
qc = I — 2sin(7r/18) for the honeycomb lattice {6,3} [9J. But generally, the substitution 
method gives us an interval in which the percolation threshold can exist. 

B. Square 

Let us now turn our attention to the star-square case shown in Fig. Mjo). By enumerating 
all the possible 345 up-sets, we find 53 different inequalities. Many of them are redundant, 
however, and we need to consider only the following inequalities: 

4p — Qp^ + Ap^ — p'^ < Aql — Aql + q^, 

with < J9 < 1 in order to find a lower bound of Pc for given qc- The first inequality is 
for an up-set {ABCD, A\BCD, B\ACD, C\ABD, D\ABC, AC\BD, A\C\BD, B\D\AC}, 
generated by {A\C\BD, B\D\AC}. The second inequality is for {ABCD, AB\CD, AD\BC} 
generated by {AB\CD, AD\BC}, and the third is for an up-set generated by {A\B\CD, 
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FIG. 2. (Color online) Solutions of the inequalities in the (a) star-triangle, (b) star-square, and 
(c) star-pentagon transformations, respectively. In (a), the horizontal dashed line represents pc = 
2sin(7r/18) ~ 0.347 296 and the vertical dashed line represents Qc = 1 — 2sin(7r/18) ~ 0.652 704. 
In (b), the dashed curve means pc = 1 — \/l — Qc to check the square lattice {4,4} (see text). 

B\C\AD, C\D\AB, A\D\BC}, respectively. As above, reverting all the three inequalities 
will yield an upper bound for pc, and the results are plotted in Fig. EJ^b). 

It is interesting to consider the square lattice {4, 4} since the star-square transformation 
transforms a square lattice with double bonds to another square lattice, rotated by angle 7r/4 
from the original lattice. Then, pc and qc should be related by pc = 1 — y/1 — qc- Our upper 
and lower bounds include this relationship over the whole region of qc- If this happened 
only within a limited region of qc, we could obtain nontrivial bounds for pc directly from 
this plot. Although this is not the case, this example shows that the substitution method 
indeed yields correct results. 



C. Pentagon 

For the star-pentagon case, the number of possible up-sets is 161 166, from which 1237 
different inequalities are found. Once again, most of them are redundant, and the set of 
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FIG. 3. (Color online) Lattice structures depicted on the Poincare disk [20]. (a) The star-pentagon 
transformation of the order-5 pentagonal tiling {5, 5} (black) leads to the order-5 square tiling 
{4,5} (gray), (b) The star-square transformation of {4,5} (black) leads to the order-5-4 rhombille 
tiling (gray), (c) The star-pentagon transformation relates the order-4 pentagonal tiling {5,4} 
(black) to the order-5-4 rhombille tiling (gray). If shifted by one lattice spacing, the order-5-4 
rhombille tiling in (c) looks the same as that in (b). 

inequalities to solve turns out to be 

5p — lOp^ + lOp^ — 5p^ + < hq^ — 5q^ + q^, 
5p^ — + p^ < q^, 
— 5p^ +p^ < hql — hql + , 

with < p < 1 when we are to find a lower bound. The first inequality is for an up-set 
generated by {AB\C\D\E, AE\B\C\D, BC\A\D\E, CD\A\B\E, DE\A\B\C}. The sec- 
ond inequality is for an up-set generated by {AB\CDE, AE\BCD, BC\ADE, CD\ABE, 
DE\ABC}. Finally, the third inequality is for an up-set generated by {A\B\CDE, 
A\E\BCD, B\C\ADE, C\D\ABE, D\E\ABC} and essentially the same as p < q^. Finding 
an upper bound is also straightforward. The results are shown in Fig. Ws). 



III. RESULTS 



A. Order-5 square tiling 

By applying the star-pentagon transformation to the order-5 pentagonal tiling {5, 5} 
with double bonds, we find the order-5 square tiling {4,5} [Fig. |3]^a)]. Note that we need 
double bonds in order to distribute five bonds to every pentagonal face. Recall that the 
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threshold pd of {5,5} is bounded as 1/4 < pH'^^ < 0.381 296 [16], which automatically 
implies 0.618 704 < pg'^^ < 3/4 by self-duality since the duality implies 

Pel +Pc2 =1, 

{m,n} {n,m} _ . 
Pc2 "T Pel — 



for lattices represented by {m,n} and {n,m} pjj. If every neighboring pair of vertices 
in {5, 5} are connected by double bonds with occupation probability p', the corresponding 
bounds of the critical threshold are located by the simple relation p'^ = 1 — v^T^^p^ as 

0.133 975 < p'll'^^ < 0.213 423, 
0.382 508 < p'll'^^ < 0.5. 

Then, solving the inequalities for the star-pentagon case, we obtain bounds for {4,5}. The 
detailed procedure is given as follows: Suppose that p'H'^^ = 0.133 975. The corresponding 
bounds are 0.133 975 < pJi'^^ < 0.413 131, whereas if p'^'^^ = 0.213 423, the bounds are 
0.213 423 < pI^'^^ < 0.527 957. Therefore, the resuhing bounds should be 0.133 975 < 
pit'^^ < 0.527 957 in total, and the same reasoning yields 0.382 508 < p^'^^ < 0.807 697. 
However, Eq. ([1]) further constrains plt'^^ as larger than or equal to 1/4. Likewise, we see 
that pH'^^ > 1/3, which implies pi2'^^ < 2/3 from the duality [Eq. ([2])]. We thus conclude 
that 

1/4 < plt'^^ < 0.527 957, 
0.382 508 < pg'^^ < 2/3. 



B. Order-5-4 rhombille tiling 

The same procedure can be repeated on the order-5 square tiling {4, 5}. The star-square 
transformation changes it to the order-5 rhombille tiling, whose face configuration can be 
denoted by V^4.5.4.5 [Fig. |3]^b)]. The face configuration means the numbers of faces at each 
of vertices around a face. The computation is similar to the above one: By putting double 
bonds between every pair of vertices in {4,5}, we see that the critical probabilities are 
bounded as 

0.133 975 < p'ii'^^ < 0.312 946, 
0.214 194 < < 0.42265. 
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When the star-square transformation is apphed, it is straightforward to obtain 

0.178 197 < p^i^-^-^-^ < 0.656 963, 
0.275 768 < p^2^-^-^-^ < 0.760 854, 

but some are no better than trivial since p^^'^'^-^ > 1/4 for coordination number n < 5, 
and the dual of V"5.4.5.4, called the tetrapentagonal tiling, has pd > 1/3 with coordination 
number 4. The bounds for this order-5-4 rhombille tiling are therefore found to be 

1/4 < p^i'^-^-^ < 0.656 963, 
0.275 768 < p^2^-^-^-^ < 2/3. 

C. Order-4 pentagonal tiling 

It is notable that the order-4 pentagonal tiling {5, 4} is also related to the order-5-4 
rhombille tiling by the star-pentagon transformation [Fig. Et^c)]. Since the duality [Eq. (|2])] 
also imposes conditions for thresholds in {5, 4} and {4, 5} as 

{4,5} {5,4} _ ^ 
Pel ^ Pc2 ~ ^ ' 

{4,5} {5,4} _ . 
Pc2 ^ Pel — ^ 5 

one could expect sharper bounds by exploiting both the relations, i.e., the star-pentagon 
transformation and the duality transformation. Unfortunately, since the star-pentagon 
transformation yields too large bounds [see Fig. Et^c)], it adds no information to the du- 
ality results, which are expressed as 

1/3 < pg'^^ < 0.617 492, 
0.472 043 < pg'^^ < 3/4, 

where 1/3 is a trivial bound from the coordination number n = 4 and p^^^^ < 3/4 is a direct 
consequence of p]^^^ > 1/4. 

IV. SUMMARY 

In summary, we have obtained analytic bounds of percolation thresholds on three hy- 
perbolic lattices by applying the substitution method to the known bounds for the order-5 
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TABLE I. Analytic bounds of bond percolation thresholds on hyperbolic lattices. 



tiling lower threshold upper threshold method 

order-5 pentagon 1/4 < p^i < 0.381 296 0.618 704 < < 3/4 Ref. [IT] 

order-5 square 1/4 < p^i < 0.527 957 0.382 508 < Pc2 < 2/3 substitution 

order-4 pentagon 1/3 < Pci < 0.617 492 0.472 043 < Pc2 < 3/4 duality 

order-5-4 rhombille 1/4 < pd < 0.656 963 0.275 768 < Pc2 < 2/3 substitution 



pentagonal tiling {5,5}. Our results are summarized in Table [H The obtained bounds are 
admittedly too broad to be very informative. But our approach illustrates how analytic 
bounds for one lattice can be made useful in estimating those for other lattices tiling a 
hyperbolic plane. Precise knowledge of Pc2 is still greatly needed in studies of percolation 
on hyperbolic lattices in general, and we hope that this approach can make further progress 
in future studies. 
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